Enhanced Elasticity and Soft Glassy Rheology of a Smectic in a Random Porous 
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We report studies of the frequency dependent shear modulus, G*(ui) = G (uf) + iG"(uj), of the 
liquid crystal octylcyanobiphenyl (8CB) confined in a colloidal aerosil gel. With the onset of smectic 
order, G' grows approximately linearly with decreasing temperature, reaching values that exceed by 
more than three orders of magnitude the values for pure 8CB. The modulus at low temperatures 
possesses a power-law component, G*(w) ~ U) a , with exponent a that approaches zero with increas- 
ing gel density. The amplitude of G' and its variation with temperature and gel density indicate 
that the low temperature response is dominated by a dense population of defects in the smectic. In 
contrast, when the 8CB is isotropic or nematic, the modulus is controlled by the elastic behavior of 
the colloidal gel. 
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Impurities and other forms of quenched disorder can 
profoundly affect the mechanical properties of materi- 
als. For instance, introducing nanometer-scale inclusions 
into polycrystalline solids is a well-developed approach 
for enhancing strength by pinning defects whose mo- 
bility dictates mechanical compliance 0. Many com- 
plex fluid phases can similarly possess topological defects 
that influence their viscoelasticity. Important examples 
are lamellar systems - such as smectic liquid crystals, 
diblock copolymers, and surfactant mesophases - that 
are characterized by one-dimensional ordering from the 
stacking of liquid-like planes. Unlike the case of crystals 
where impurities can be frozen into the lattice, intro- 
ducing quenched disorder into such fluids is a technical 
challenge. In some cases, annealed disorder in the form 
of colloids suspended in lamellar systems can anchor line 
defects into a crosslinked network that enhances the elas- 
tic response 0, • 

One systematic approach for introducing quenched dis- 
order to a smectic is to confine the fluid in a random 

porous medium HHE0,SH[Ii[l3- Study ° f Such 
systems has largely concerned the influence of the en- 
vironment on the phase behavior and smectic ordering. 
Theory and experiment have concluded that the smectic 
phase is highly sensitive to the quenched disorder. An 
interesting example is a smectic confined by an aerosil 
gel, which is a compliant hydrogen-bonded gel formed 
from nanometer-scale silica colloids [jj la 0, la, m In- 
consistent with theory 0, > x-ray scattering on liq- 
uid crystal - aerosil composites has revealed that the gel 
introduces random field effects that replace the thermo- 
dynamically sharp nematic to smectic transition with the 
growth of short-range smectic correlations 0,0,01 ■ The- 



ory further predicts that quenched disorder impacts the 
mechanical behavior of smectics by introducing enhanced 
anomalous elasticity |ll|. We present rheometry stud- 
ies that characterize how the viscoelastic properties of a 
smectic are influenced by such disorder. Using weakly 
bonded aerosil gels enables us to isolate the mechanical 
response of the smectic. We find that the quenched dis- 
order greatly increases the shear modulus over that of 
the pure smectic and leads to weak power-law behavior 
characteristic of soft glassy rheology. 

Samples were prepared by forming aerosil gels within 
the liquid crystal octylcyanobiphenyl (8CB) following es- 
tablished procedures J5|- The 8CB (Frinton Labs) had a 
quoted purity of 99.4%. Pure 8CB undergoes an isotropic 
to nematic transition at Tffi — 313.98 K and a nematic 
to smectic-A transition in which the molecules form a 
density wave with period do = 3.4 nm at Tna = 306.97 
K |5|. The aerosil (DeGussa Corp., type 300) consists of 7 
nm Si02 particles that are hydrophilic due to a high den- 
sity of surface hydroxyl groups. Degassed 8CB and dried 
aerosil were mixed with high purity acetone and soni- 
cated for several hours to achieve a uniform dispersion. 
The mixtures were then heated to 320 K to evaporate 
the acetone slowly. During this process the aerosil aggre- 
gated via hydrogen bonding. The resulting composites 
were uniform, thixotropic solids comprised of 8CB per- 
meated by a colloidal gel. Gels formed in this way have a 
fractal dimension near 2.1 and a broad distribution 
of void sizes with a mean pore chord Iq that varies with 
aerosil concentration as lo = 6.7 j ps nm, where ps is the 
aerosil density in grams Si02 / cm 3 8CB 8] . Samples were 
prepared with ps — 0.035 to ps = 0.10, corresponding 
to Iq ~ 56do to Iq w 19g?o arid aerosil volume fractions 
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of 0.016 to 0.043. For ps > 0.10 the composites change 
from a pasty to a brittle consistency that makes rhcom- 
etry measurements difficult. 

The gel affects both the liquid crystal's nematic and 
smectic behavior. High porosity gels weakly perturb the 
orientational ordering, as evidenced by nematic correla- 
tion lengths that exceed Iq by a factor of 30 or more for 
ps < 0.10 fiui(. The observed random field behavior of 
the smectic Q is understood by considering how the ran- 
domly positioned gel strands impose steric hinderance on 
the liquid crystal molecules, effectively forcing the phase 
of the smectic mass density wave to different values at 
different points in space. The sensitivity of the smectic 
ordering to the gel is reflected in the smectic correlation 
lengths, which are nearly an order of magnitude smaller 
than the nematic correlation lengths [7J. 

The frequency-dependent shear modulus, G*(uS) = 
G'(id)+iG"(uj), of the 8CB-aerosil composites was deter- 
mined using a stress controlled rheometer (Paar Physica 
MCR300) in a cone and plate geometry equipped with 
a Peltier device that maintained temperature stability 
to within ± 0.01 K. Temperature gradients across the 
sample were estimated at less than 0.2 K. Figure 1(a) 
shows the storage modulus G'(lo = 1 rad/s) as a func- 
tion of strain amplitude 70 for ps = 0.10 at temperatures 
above and below the smectic transition. When the 8CB 
possesses smectic order, the response depends weakly on 
70 to the smallest strains. However, measurements of 
G*(lo) with different 70 showed only this slight variation 
in magnitude; all salient features, particularly the depen- 
dence on u) and temperature, were independent of 70 for 
7o < 0.007. The results presented below were obtained 
with 70 = 0.003. 

Figures 1(b) and 1(c) show G'(lo) and the loss modulus 
G"(co) respectively at several temperatures for ps = 0.10. 
Above 307 K the modulus obeys trends of a weak elastic 
solid, with G'(lo) constant and G'(lo) » G"(ui), and dis- 
plays little temperature dependence. In particular, the 
isotropic to nematic transition near 314 K has no measur- 
able effect on G*(ui). We associate this high temperature 
mechanical response with that of the aerosil gel in a low 
viscosity liquid crystal solvent. Measurements of aerosil 
gels in the simple organic liquid dibutyl phthalate give 
quantitatively similar G* (lo) , supporting this interpreta- 
tion. The storage modulus G g above the smectic tran- 
sition temperature varies with aerosil density roughly as 
G g ~ j o|' 8±0 ' 8 , consistent with behavior expected for a 
colloidal gel with fractal dimension near 2.1 [l^. With 
the onset of smectic ordering in the 8CB at low temper- 
atures, G' (lo) increases dramatically and obtains a clear 
frequency dependence. We attribute the growth in G*(u>) 
to contributions from smectic defects. The solid lines in 
Fig. 1(b) are the results of fits to a power-law contribu- 
tion plus a frequency-independent background 



G'(u>) = G g + G s + (i s io c 



(1) 



where G g is the contribution from the gel, assumed 
independent of lo and temperature, and G S (T) is the 
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FIG. 1: (a) Storage modulus G' of 8CB confined in an aerosil 
gel with ps = 0.10 g/cm 3 at lo — 1 rad/s as a function of strain 
amplitude 70 below (solid circles) and above (+) the smectic 
transition temperature, Tna- (b) G' '(lo) at several temper- 
atures for 70 = 0.003. Lines are results of fits to Eq. (1). 
(c) Loss modulus G" (lo) below (open circles) and above (Xs) 
Tna- The solid circles display the difference between the two 
spectra, and the line is the result of a power-law fit. 



frequency-independent contribution from the smectic 
phase. This form describes accurately the storage mod- 
ulus at all densities of aerosil and at all temperatures. 

The inset to Fig. 2 displays the amplitudes for the 
power-law and frequency-independent contributions, f3 s 
and G s , as a function of temperature for ps — 0.10. These 
parameters show qualitatively the same temperature de- 
pendence for all ps except that the ratio fl s /G s increases 
with increasing ps- The exponent a characterizing the 
power-law dependence displays no systematic tempera- 
ture dependence below the pseudocritical point marking 
the onset of static short-range smectic order 0, HJ but 
does vary systematically with ps- Figure 2 displays the 
value of a averaged over temperature as a function of 
ps wherein a decreases from a w 0.25 at ps — 0.035 to 
a w 0.1 at ps = 0.10. Such weak power-law behavior 
is characteristic of a variety of complex fluids including 
foams, emulsions, particulate suspensions and slurries. 
In a recent model Sollich and coworkers jig provide a 
unifying theoretical framework for this "soft glassy rheol- 
ogy" , arguing that the response is a general consequence 
of structural disorder and metastability. In this model 
a serves as an effective noise temperature, with systems 
enterin g a glass transition as a — ► 0. (The exponent x 
in Ref. [lq is x = a + 1.) Thus, the variation of a with 
ps points to increasingly glassy dynamics in the smectic 
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FIG. 2: The exponent a characterizing the power-law con- 
tribution to G'(cu) (solid circles) and G"(a>) (open circles) 
for smectic 8CB confined by an aerosil gel as a function of 
gel density. The inset displays results for G s (open circles) 
and I3 S (filled circles) corresponding to the amplitudes of the 
frequency-independent and power-law contributions to G'(u>), 
respectively, for 8CB in a gel with ps = 0.10 g/cm 3 at tem- 
peratures below the onset of smectic order in the 8CB. 



with increased quenched disorder. This notion of glassy 
behavior further connects the rheology with the slow dy- 
namics observed with dynamic light scattering for 8CB in 
aerosil 17]. However, we note the simultaneous growth of 
the frequency-independent component G s implies some 
portion of the 8CB is frozen for all p$. 

As expected for soft glassy rheology the weak 

power-law behavior also appears in the loss modulus. 
Figure 1(c) displays G"(lo) at low temperature in the 
smectic phase of 8CB and at high temperature in the 
nematic phase, where contributions from the gel pre- 
sumably dominate. With the high temperature modu- 
lus subtracted from the low temperature spectrum, the 
remaining smectic contribution, shown by the solid sym- 
bols in Fig. 1(c), displays power-law behavior. As shown 
in Fig. 2, fits to such results, G"(uj) ~ G'(u>) ~ uo a , pro- 
vide values for a that agree closely with those obtained 
from G'(uj). 

The amplitude and frequency dependence of G'(u>) for 
smectic 8CB in aerosil differ from those of other un- 
aligned lamellar complex fluids. A wide range of systems 
including lyotropic and thermotropic smectics (without 
quenched disorder), surfactant-based phases, and diblock 
copolymers exhibit a frequency response with two con- 
tributions: G'(uj) — Go + fid^ 1 / 2 . The plateau modu- 
lus Go arises from the elasticity of static defects while 
the frequency-dependent term results from the bulk re- 
sponse of t he ma terial. Specifically, ftdUJ 1 / 2 with (3d = 
(71-/24^)^(57?), where B is the compression modulus 
and r] is an effective viscosity, describes the response of 
lamellar regions with the layer normal oriented so that 
strain necessitates layer compression . While the rel- 
ative strength of these terms varies with material param- 
eters and defect structures, most studies on unaligned 
lamellar systems capture this characteristic response in 
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FIG. 3: (a) Storage modulus G'(T) and (b) loss modulus 
G"(T) at CJ = 1 rad/s for 8CB in an aerosil gel with ps = 
0.035 g/cm 3 (triangles), 0.05 g/cm 3 (squares), 0.07 g/cm 3 
(diamonds), and 0.10 g/cm 3 (circles), (c) The ratio G"/G' 
scaled by the exponent a for the four densities. 



their experimental frequency range. In particular, exper- 
iments on pure unaligned 8CB display a plateau modu- 
lus Go ~ 100 Pa and a ui 1 ^ 2 contribution with (3d ~ 30 
Pa-s 1 / 2 [l9l I20I]. consistent with the room temperature 
values B w 2 • 10 6 Pa and 77 w 0.1 Pa-s. The response 
we observe for 8CB in aerosil is more than three orders of 
magnitude largerth&n G'(u>) of pure 8CB. The greatly en- 
hanced modulus in the presence of disorder implies that 
contributions from the bulk fluid are negligible. Instead, 
we infer that the shear response derives from a dense 
population of defects that form in the smectic with the 
destruction of smectic correlations by the disorder. Con- 
straints imposed on defect motion by the confining gel 
provide rigidity that dominates the glassy response to 
shear. 

Evidence for soft glassy rheology and the approach to a 
glass transition with increasing ps comes not only from 
the weak power-law behavior of the modulus but also 
from the ratio G"/G'. Figures 3(a) and 3(b) show the 
temperature dependence of G' and G", respectively, at 
u = 1 rad/s for varying p.g. Within the soft glassy rhe- 
ology model G"/G' ~a0. Figure 3(c) shows G"/G' 
scaled by a. At temperatures below the pseudocritical 
region, the ratios collapse onto a single scaling function. 
This scaling demonstrates how the single parameter a 
accounts for the variations in mechanical response with 
changing ps- Thus, a direct correspondence exists be- 
tween the strength of random fields in the system, as 
parameterized by ps 0- Q ■ and the exponent a that is 
central to soft glassy rheology. 

The variation of G' with ps indicates that smectic line 
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FIG. 4: Smectic contribution to G' near the pseudotransition 
temperature that marks the onset of smectic order normalized 
by its value at the reduced temperature e = 5 • 10 for 8CB in 
an aerosil gel. Included are G'(w = 1) for ps = 0.035 g/cm 3 
(triangles), 0.05 g/cm 3 (squares), 0.07 g/cm 3 (diamonds), and 
0.10 g/cm 3 (circles) as well as G'(u> = 0.1) (crosses) and 
G'(u) = 40) (X's) for 0.07 g/cm 3 . The solid line represents a 
linear dependence on e. The inset displays the values of G s at 
292.6 K with the solid line showing the fit result G s ~ Ps' 9±0 ' 2 . 

defects dominate this soft glassy rheology. In analogy 
with rubber elasticity, the elastic response of a static de- 
fect network should vary as G s » t/cP where r is the 
defect line tension and d is a typical spacing between de- 
fects 0- A plausible scale for d is set by the porosity of 
the gel. Assuming d f=a £ ~ Iq, where £ is the smectic cor- 
relation length and Iq is the gel's mean pore chord p|0, 
we expect G s ~ ^ 2 ~ p|. The inset to Fig. 4 displays G s 
at 292.6 K, well below the nematic to smectic transition, 
as a function of ps. The solid line in the inset represents 
the result of the best fit to a power law, G s ~ Ps 9±0 ' 2 , 
consistent with this expectation. 

The temperature dependence of G' similarly indicates 
defects dominate the mechanical response. Figure 4 
shows the smectic contribution to G' as a function of 
reduced temperature, e = (Tna — T)/Tna, normalized 



by its value ate = 5-10~ 3 [2l|. Included in the plot are 
data from Fig. 3(a) covering G'(u> — 1, e) for four aerosil 
densities as well as G'(u> = 0.1, e) and G"(w = 40, e) tra- 
ps = 0.07. The data collapse onto the solid line which 
displays the relation (G'(e) — G g ) ~ e z with z = 1. This 
scaling demonstrates that the smectic contribution to the 
mechanical response is associated with critical phenom- 
ena of the nematic to smectic transition. As mentioned 
above, analogy with rubber elasticity predicts G' ~ r, 
where r is the defect line tension. For screw dislocations, 
which are prominent in smectics j2^, t = Bb 4 /128Tr 3 r 2 , 
where b — mdo is the Burger's vector of integer strength 
to and r c is the defect core radius |2^|. For 8CB near 
T NA , B ~ e 0A H| and r~ 2 ~ ^ 2 ~ e 5 , where V 
is the smectic order parameter 0. Thus, r ~ e 9 in 
close agreement with zwlwe observe for G", suggesting 
that screw dislocations play a predominant role in dictat- 
ing the shear response. A more quantitative assessment 
of the prominence of screw dislocations depends on the 
typical Burger's vector strength. For a moderate value, 
to = 4, along with the low temperature values B w 10 7 
Pa [24(, r c « 1 nm, and G s from the inset of Fig. 4, the 
relation G s w r/d 2 leads to a plausible defect spacing 
d « 0.5/o, independent of aerosil density. 

Regardless of the precise connection between the tem- 
perature dependence of G' and quenched disorder, the 
observed enhanced elasticity and soft glassy rheology of 
8CB confined in aerosil provides a new perspective on 
effects of quenched disorder in complex fluids. Unlike for 
many soft glassy systems, the defects that drive this rhe- 
ology provide a clear physical origin for the disorder and 
metastability in the smectic. Furthermore, the quenched 
disorder that the gel imposes on the smectic has been con- 
vincingly characterized in terms of random fields . 
Thus, this mechanical behavior connects a classic model 
of quenched disorder, from random fields, with a new 
paradigm of glassy behavior in complex fluids. 
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